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CHARACTERISTICS OF TUNNELING P-N JUNCTIONS 


By John B. Hopkins 
Electronics Research Center 


SUMMARY 


Procedures have been developed for the calculation of current -voltage 
characteristics of tunneling p-n junctions in semiconductors, given material 
constants and impurity -doping parameters. Accuracy higher than that 
acquired in similar efforts has been attained. When a digital computer is 
used, numerical differentiation permits calculation of parameters relevant to 
microwave mixing and detecting applications for tunnel and backward diodes 
in any elemental or compound semiconductor. Typical results are indicated. 


INTRODUCTION 

This report concerns research done in the fabrication and characterization 
of backward diodes made from a variety of III-V semiconductors. They are in- 
tended for use in microwave mixing and detecting applications. The research 
effort requires estimation of tunnel diode characteristics in a number of mate- 
rials and for varying impurity concentrations. For example, it is useful to be 
able to estimate figures of merit which involve the first and second derivatives 
of the i-v curve, and these can be meaningful only with a fairly accurate know- 
ledge of i(v). Determination of the theoretical i-v characteristics and other 
parameters of a p-n junction tunnel diode is a problem not amenable to an ana- 
lytical approach. 

A number of authors (refs. 1, 2, 3 and 4) have used greatly simplified 
models to obtain various numerical or analytical results, but--in addition to 
their very approximate nature- -the se solutions are limited in the information 
provided. This report describes results obtainable by use of an analysis sig- 
nificantly more detailed than that conventionally used. These calculations, 
facilitated by use of a digital computer, can be made to a significantly higher 
degree of accuracy than those previously reported. 

It should be noted, however, that the calculations, though detailed, are 
ultimately based on the conventional physical model, and are therefore sub- 
ject to errors resulting from the inadequacies of that representation. These 
weaknesses will be indicated at relevant points in the detailed description of 
the model and the calculations. 



THE PHYSICAL MODEL 


The model used is a conventional p-n junction. Perfect crystal structure 
is assumed, with an impurity concentration which changes abruptly from N^, 
at energy -E^ (n- region), to N a , at energy E a (p- region). It should be noted 
that no real junction can be truly abrupt, and this assumption may well lead to 
inaccuracies when applied to tunneling junctions. For these, the "abruptness” 
condition (that the depletion layer be much wider than the region of changing 
impurity concentration) is seriously violated. Nor are defects of any type con- 
sidered, so that this program does not show excess current. While this is a 
definite limitation for application to tunnel diodes (the peak- to- valley current 
ratio cannot be estimated, for example), it is forced by the theoretical com- 
plexity of the problem. 

Further, it is of less relevance to backward diode studies, where opera- 
tion is generally at very low bias, and excess current is small. The band 
structure is taken to be particularly simple- -parabolic, andisotropic in k. 

Also, only direct tunneling transitions are considered. This could be a seri- 
ous restriction, since both germanium and silicon are indirect-gap materials. 
However, the principal error will be only in current densities, so much useful 
information can be obtained for these materials as well. (Most III-V semi- 
conductors of current interest are essentially direct-gap and accurate esti- 
mates can be obtained for them. ) 

Tunneling current is a very sensitive function of tunneling distance, and 
considerable care is taken in determining this parameter. The various mate- 
rial constants are specified for the temperature assumed, and are taken to be 
independent of applied voltage and Fermi levels, and as constant throughout 
the material. However, variation of the diffusion constant with the position 
of the Fermi level is included. 

Discussion of tunnel diodes necessarily requires consideration of degen- 
erate and near- degenerate conditions, and difficulties do arise--for example, 
in determination of the Fermi levels. Fortunately, some uncertainty in this 
does not seriously affect results in the region of principal interest. Further, 
an approximation has been developed which seems sufficiently accurate. 

The required data for the material includes band gap, dielectric constant, 
lattice constant, and hole and electron mobilities, effective masses, and re- 
combination times. The diffusion constants are calculated from the mobilities, 
using the Einstein relation (D = kTp/q), corrected for degeneracy. The diffu- 
sion lengths are determined from the diffusion constants and recombination 
times 


[ i = 7T ]. 
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In addition, the impurity doping parameters are read in. These include both 
concentrations and energy levels. Variation of mobility caused by increased 
impurity scattering for increasing doping density is incorporated by use of the 
Dingle- Brooks expression for impurity mobility. 

As the program is intended to provide information about phenomena asso- 
ciated with low voltages only, no breakdown mechanisms have been included in 
the diode model. 
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DETAILS OF THE CALCULATIONS 


The first step in determining the electrical characteristics of a semi- 
conductor is calculation of the position of the Fermi level. For a p-n diode, 
this is done separately for the p and n regions. The basic equation used is 
that of charge neutrality: 


electrons in 1 

holes in 1 

’ ionized"! 

’un -ionized 1 

.conduction bandj 

_ valance band J 

donors J 

[acceptors J 


= 0 ( 1 ) 


The assumption is made that in the n region N = 0, and in the p region N<J = 0. 
(These restrictions could easily be removed, but have been made to avoid un- 
necessary complication of the program. ) Through entirely conventional man- 
ipulations, equation (1) may be written: 
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An analogous equation describes the p region. All three terms in equation (2) 
are functions of Ef, and no closed analytic expression can be obtained. A 
simple iterative procedure is used (with a separate subroutine for the com- 
puter calculation) to determine the Fermi levels. For E less than -5, the 
Ehrenberg approximation 

^F ermi integral = 2 tt^^ 



is used; otherwise numerical integration is carried out with increments of 0. 1 
over the range from zero to 


25 + 


6 f 
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(A change of variable u = e is made to cause the integrand to vanish more 
quickly. ) The iteration for E, consists of successive halving of an interval 
originally chosen to include the expected value. Any desired accuracy can 
be obtained; the value used for the uncertainty in Ef is 0. 01 (. 00026 eV). 


When calculations of Ef for the n and p regions have been made, carrier 
concentrations are immediately determined from 
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If the material mobility is known, conductivity may then be calculated as 

l h * 


^ = n qp e + pqp^ 


Degeneracy (impurity concentrations such that the Fermi levels are 
within the n- region conduction band and the p- region valence band) presents 
some difficulty. At these impurity concentrations, the assumption that all 
impurities are at a particular energy level, with no interaction, is signifi- 
cantly in error. The theoretical situation is quite involved, and no simple, 
generally valid result is available. However, at this point the dependence of 
device characteristics becomes somewhat less critical. An emprical ap- 
proximation -- f(Ef) = 1/3 for E£ above the donor level or below the accep- 
tor level -- seems to agree fairly well with relevant published data (e.g. , 
ref. 4). (This gives a continuous function, since f (E v ) = f n (E c ) = 1/3.) 
Such an approximation is unlikely to be accurate at strong degeneracy, but 
that is a situation rarely found. For backward diodes (the principal subject 
of this study) only slight degeneracy occurs. Mobility is affected markedly 
by high impurity concentrations, and a correction is necessary. The one 
chosen is that of Dingle and Brooks (refs. 5, 6), for which the mobility \i ^ 
caused by impurity scattering is given by 
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with a the mobility at low doping, caused principally by lattice scattering. 
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One further correction for degenerate circumstances involves the diffu- 
sion constants, which are also somewhat affected by heavy doping. Again a 
small empirical correction is made: D = (1 + 2E^) D q within the conduction 

and valence bands, where D Q is the value at low doping. 

The next step is computation of device behavior under applied bias. All 
voltage drop is assumed to take place across the junction, and -- following 
generation of the specified applied voltages -- the thickness and shape of the 
depletion region are determined. Again, this is a conventional calculation, 
based on the assumption of an abrupt junction, and is outlined in Appendix I. 
Junction capacitance (per unit area) is then calculated from the depletion 
layer width and the dielectric constant of the material. For a specified junc- 
tion area, junction (barrier) capacitance C-g is then known, as is the spread- 
ing resistance, Rg, which is determined from the semiconductor conductivity: 
Rg = 1/4 err, where r is the junction radius. 

As indicated in Appendix I, the actual tunneling distance is different 
from the junction width, though related to it, and this is used in the next 
section where tunnel current is calculated. 


The tunnel current across a p-n junction is generally given by 
E 

I t = J p c (E)p v (E)f n (E)[l -f r (E)]P t dE, E v >E c (3) 
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The tunneling probability is calculated from quantum mechanics, with a 
triangular energy barrier assumed, and is usually taken as 
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where F, the average field across the junction, is taken as 
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with the tunneling distance t = t(E), as given in Appendix I. 

The numerical integration (according to Simpson's rule) is carried out 
by dividing the energy interval into 20 parts, which is found to give results 
not significantly different from division into 100 parts. As indicated before, 
the energy bands are assumed to be parabolic, so 
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and 


p c (E) a (E - E ^) 2 
P V (E) ct (E v - E ) 2 


within the bands and zero elsewhere. 

For reference, the normal diode current is also calculated. This is 
given by 


I = I s [exp(qV ap /kT) - 1] (5) 

with I determined from carrier concentrations, and diffusion constants and 

g 7 

lengths: 


I = q[D n /l -F D p / L. ] . 

s ^ n o n p *o p 

p r n 

The validity of this result is somewhat limited, as diffusion lengths are 
seldom known with any precision, and, in addition, other mechanisms often 
add considerably to this current. However, it is generally not very signifi- 
cant in the region of interest (low applied voltages), particularly when sub- 
stantial tunneling occurs. 

The peak voltage is determined by a simple iterative scheme. Tunnel 
current is calculated for increasing values of voltage until two successive 
currents show a decrease. The procedure then repeats from the last 
increasing point, with a reduced voltage increment. After several such 
cycles, an accuracy of better than 0. 1 mV can be obtained. Thus, for each 
specified value of applied voltage, the program calculates resulting tunnel 
current, diffusion current, junction width and capacitance, and "average" 
tunneling distance (see Appendix I). The value of current associated with 
peak voltage is then used to normalize the tunnel current values. In addition, 
useful information concerning device properties can be obtained from the 
first and second derivatives of the current- voltage characteristic. For 
example, current sensitivity |3 0 is given by f M /f ! where i = f(v). Numerical 
differentiation is difficult to accomplish with significant accuracy, particu- 
larly for higher derivatives. The method used here involves calculation of 
current- voltage data points displaced from the point of interest by plus-and- 
minus 1- and 2- volt increments. These are used (as shown in Appendix II) 
to approximate f r (v) and f rt (v). The size of the voltage increment to be used 
must be a compromise, as the natural choice of a small increment results 
in a set of data points so close together that differences between them become 
very small, and even the eight- place accuracy of a computer is inadequate. 
The optimum size has been determined by comparison of different trials, 
taking as a criterion that variation should be smooth from point to point. 

An increment of 1 mV has proven satisfactory and should prove compatible 
with experimental techniques. 


The physical model used does not permit determination of actual cur- 
rents or current densities with meaningful accuracy. However, published 
experimental data can be used to obtain approximate calibration constants 
to relate the calculated values to actual currents. With these estimates, 
the barrier resistance R-^(l/f I ) is easily calculated, and -- with the capaci- 
tance Cg and spreading resistance Rg obtained previously - -a number of 
figures of merit useful for microwave device characterization can be found: 


P = tf'/f 1 

[1+ (R s /R b )]_ 1/2 
c - 2 ttC b (R s R b ) 1/2 


< 1 +R s /R B> [ 1 + ( f o /f c^ 2] 


M 
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Running time on a digital computer for these calculations can be quite 
short. On a Univac 1108, compilation requires 4 seconds, and the deter- 
mination of Fermi levels, carrier concentrations, peak voltage and current, 
and so forth, takes approximately 0. 5 second. Several hundred i-v data 
points may be calculated in 1 second. 
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SUMMARY AND CONCLUSIONS 

The calculations described above may be thought of as comprising two 
major parts -- computation of Fermi levels, and determination of resulting 
material and diode parameters . The accuracy of the Fermi level calcula- 
tion is difficult to check, but agrees fairly well with published estimates. 

One check is comparison with the data compiled by Irvin (ref. 7) showing 
resistivity versus doping for n-type silicon. His values represent an aver- 
aging of published experimental data, with 90 percent of the actual points 
falling within 10 percent of his curve (below = 10*9). The match of 
calculated resistivity with Irvin* s curve is shown in Figure 1, and is quite 
good up to 5 x 10^9 9 at which point impurity scattering and conduction pro- 
duce a situation too complex to treat usefully in any simple fashion. 

Typical calculations of Fermi level as a function of impurity concentra- 
tion, near the band edge, are shown in Figure Z. The current- voltage 
characteristics for a particular diode is graphed in Figure 3 where two 
calculated curves are compared to measured values for a GaSb unit. One 
calculation is simply 

E^ 

J P C (E) p v <E)f n (E) [l-f p (E)]dE 
E 

c 

as used by Esaki (ref. 4) and others, and for which a good numerical ap- 
proximation has been developed by Bates (ref. 8). The other computation 
includes the tunneling distance and probability as variables within the 
integral. The latter approach is seen to offer substantially higher accuracy, 
matching quite well until excess current, probably associated with tunneling 
between impurity states, becomes significant. (N^ and N a were chosen in 
each case to match the peak voltage of the experimental curve. ) 

Figures 4 and 5 are typical computer printouts. Figure 4 illustrates 
simple i-v curve data, while Figure 5 is a sample of the figure of merit 
calculations . 

The calculations that have been described are not presented as an 
exact solution of high accuracy. Indeed, many approximations and weak- 
nesses have been indicated. However, computations based on this model 
have been found to be of use, providing significant improvement over the 
simpler approaches generally taken in articles and textbooks to obtain es- 
timates of both numerical values and the relation between different vari- 
ables. For example, the improvement in determination of current- voltage 
characteristics should make possible meaningful comparison of various 
semiconductors in terms of their suitability for fabrication of backward 
diodes to be used for microwave detection and mixing. Further, useful 
information can be obtained concerning optimum doping and performance 
as a function of temperature. 
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Degeneracy (mV) 






Normalized Current 


In summary, this formulation permits rapid and easy estimation of ap 
proximate device characteristics, with significantly better accuracy than 
has been possible with previous simpler approaches. 
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Applied Voltage (mV) 

Figure 3. - -Comparison of calculated and experimental tunnel diode 
characteristics (GaSb) (matched for peak current and voltage; 
equal degeneracy for both bands) 


germanium 


E(G) = 

• 07 

ME) = .22 

M(H)= .30 

MU(N>= 3800.00 

MU(P)= 1800.00 

D(N> = 

98.22 

D(P)= 46,53 

L(N)= 70.33 

L(P)= 41,43 

T(N)= .50-03 T(P)= .50-03 

Dielectric constant: 

= 16.00 

LATTICE CONSTANT: 

5.66 ANGSTROMS 


TEMP (K 

•) = 300.0 

JUNCTION 

AREA A. 50 X 2.50 MICRONS 

= .625-07 SQ. CM 

. MOBILITY ( IMPURITY) N= 534. P= 

DUPING 

N ( D ) = 

,600 + 20 

N(A)= .600+20 




E(D) = 

-.030 EV. 

E(A)= .030 EV. 




EF(N)= 

.117 EV. 

EF(P)= -.002 EV. 

FERMI 

FUNCTION FROZEN AT EC(N) = -.030 EC(P) = 


PtAK CURRENT = .1052-04 

PEAK VOLTAGE = 50.56 MV. 

BAND OVERLAP = 

199.7 MV. 

IC/CJ 

= .084973657 

N0(N)- 

.20+20 

PO(P) = 

.20+20 






PO(N)= 

.24+06 

NO(P) = 

.59+06 






RH(J(N1 = 

,564-03 OHM-CM. 

RHOIPU 

.869-05 OHM-CM. 






RS(N)= 

1.04 OhMS 

RS(P) = 

1.54 OHMS 






VOLTAGE 

TUN. CURRENT 

DIF, CURRENT TOT. CURRENT 

TUNNEL DISTANCE 

JUNCTION WIDTH 

JUNCTION CAPACITANCE 

-.050 

-3.9867 

-.885-17 

-.419-04 

37.6 

A. 

73.7 

A. 

.1204 

-.040 

-2,6788 

-.815-17 

-.303-04 

37,9 

A. 

73.3 

A. 

.1210 

-.030 

-1,9373 

-.711-17 

-.204-04 

38.3 

A. 

72.9 

A. 

.1217 

-.020 

-1,1518 

-.557-17 

-.121-04 

36.3 

A. 

72.5 

A. 

,1224 

-.010 

-.5101 

-.332-17 

-.536-05 

39.0 

A. 

72.1 

A. 

.1231 

-.000 

.0000 

,000 

,000 

39.4 

A. 

71.7 

A. 

.1238 

• 010 

.3905 

,489-17 

.411-05 

39.8 

A. 

71.3 

A. 

.1245 

.020 

.0735 

.121-16 

.708-05 

40.2 

A. 

70.8 

A. 

.1252 

.030 

.8611 

.227-16 

.906-05 

40.7 

A. 

70.4 

A. 

.1260 

.040 

.9657 

.363-16 

.102-04 

41.1 

A. 

70.0 

A. 

.1267 

« 05U 

.9999 

,613-Id 

.105-04 

41.6 

A. 

69.6 

A. 

.1275 

• 060 

.9767 

.951-16 

.103-04 

42.2 

A. 

69,2 

A. 

.1283 

.070 

.9091 

.145-15 

.956-05 

42.7 

A. 

68,7 

A. 

,1291 

.080 

.0104 

.218-15 

.852-05 

43.3 

A. 

68,3 

A. 

.1299 

• 09U 

.6929 

.326-15 

.729-05 

43.9 

A. 

67.9 

A. 

.1307 

.100 

.5662 

.485-15 

.598-05 

44.6 

A . 

67.4 

A. 

.1316 

.110 

.4458 

.719-15 

,469-05 

45.4 

A. 

67.0 

A. 

.1324 

.120 

,3333 

.106-14 

.351-05 

46.2 

A. 

66.5 

A. 

,1333 

.130 

.2359 

.157-14 

.248-05 

47,0 

A , 

66.1 

A. 

.1342 

.140 

.1564 

.232-14 

.164-05 

48.0 

A. 

65.6 

A. 

.1351 

.150 

.0954 

,342-14 

.100-05 

49.1 

A. 

65.2 

A. 

.1361 

.160 

.0520 

.504-14 

,547-06 

50.4 

A. 

64.7 

A. 

.1370 

.170 

.0239 

.742-14 

.252-06 

49.1 

A. 

64.3 

A. 

.1380 

.180 

.0082 

,109-13 

.862-07 

53.7 

A. 

63,8 

A. 

.1390 

.190 

.0014 

.161-13 

.147-07 

56.3 

A. 

63.4 

A. 

.1400 


Figure 4. - -Simple i-v curve data 
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Figure 5. --Figure of merit calculations 
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APPENDIX I 


DETERMINATION OF DEPLETION LAYER AND TUNNELING 

THICKNESS 


The quantities necessary for the determination of depletion layer and 
tunneling thickness are indicated in Figure A- 1 . The zero of potential is 
taken as the bottom of the n- region conduction band, far from the junction. 
The doping is assumed to change abruptly from N<j (n region) to N a (p re- 
gion) at x = 0. It is convenient to define V n (x) and Vp(x) as the potential 


v 



Figure A-l. --p-n junction showing relative energy levels and distances 

at the bottom of the conduction band in the n and p regions, respectively. In 
completely conventional fashion, Poisson's Equation is integrated twice to 
give 
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The terms dp and d n represent the x- values for which space charge (and 
therfore the electric field) is taken as zero. This condition defines the 
boundaries of the depletion region. and Vq are integration constants 
and must be equal, since Vp(0) = V n (0) = = V^. Further 
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From Figure A- 1 , it can be seen that 
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An equal amount of (space) charge must exist on either side of the junction, 
so that 
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and 



with the total width of the depletion region d = d n + dp, and the junction ca- 
pacitance C = Ae/d. 1 

The tunneling distance, for carriers of a particular energy, is x^ - Xp, 
with Xp generally negative. This distance exists only in the range 6v > 0. 

For-d^< x < 0, the expression for the top of the p- region valence band is: 

qN a 2 N d d n 

V (x) - E = (x + d x) + q — - — y - E , 

P g c P H « 2 g ’ 

and for the bottom of the n- region conduction band: 
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For 6 vq, a particular value of V, with 0 < Sv Q < 6v, x^ and Xp are implicitly 
determined from 
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or rewriting: 
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(The correct signs for the radicals are determined by consideration of lim- 
iting cases. ) Thus: 



For specified materials and doping, the degeneracy (&Vp and Sv n ) is known, 
and 6v is linearly related to the applied voltage. To calculate the tunnel cur- 
rent, one can now integrate over V from zero to &v, with tunnel distance t 
being a known function of potential. An "average" tunnel distance can also 
be determined. This quantity is somewhat arbitrary; in this program it is 
defined by 


t = 
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J * t(v) dv. 
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APPENDIX II 


NUMERICAL DETERMINATION OF FIRST AND SECOND 

DERIVATIVES 


Consider a function f(x), known at a set of equally spaced points x_^, 
x -(k- 1 )> ••• x -l» x o> x l> ••• x k» witl1 f k = f ( x k^* ° ne can define "central 
differences" of various order as follows: 
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As shown in standard references on numerical analysis (ref. II- 1), the first 
and second derivatives of f, evaluated at x Q , are then given by 
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+ terms of higher order 
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+ terms of higher order 


with h 


= x k - *k-i- 


These expressions are used in the diode computer program, with h = 1 mV. 
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